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Monopole-antimonopole pairs connected by strings and monopole- 
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. string networks with N > 2 strings attached to each monopole can be 

. formed at phase transitions in the early universe. In such hybrid defects, 

C . monopoles accelerate under the string tension and can reach ultrarelativis- 
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Qs, . accelerating monopoles. For monopoles with a chromomagnetic charge, we 

r-j 

(-^ . also discuss the high-energy hadron production through emission of virtual 



tic Lorentz factors, 7 » 1. We study the radiation of gauge quanta by 



gluons and their subsequent fragmentation into hadrons. The relevant pa- 



CZ2 . 

c3 rameter for gauge boson radiation is M/a, where M is the boson mass and 



a is the proper acceleration of the monopole. For MCo, the gauge bosons 
can be considered as massless and the typical energy of the emitted quanta is 
E ~ 7a. In the opposite limit, M 3> a, the radiation power is exponentially 
suppressed and gauge quanta are emitted with a typical energy E ~ 7M 
in a narrow range AE/E ~ (a/M) 1 / 2 . Cosmological monopole-string net- 
works can produce photons and hadrons of extremely high energies. For 
a wide range of parameters these energies can be much greater than the 
Planck scale. 
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I. INTRODUCTION 



Monopoles connected by strings can be formed in a sequence of symmetry breaking 
phase transitions in the early universe [T]Q]. A typical sequence of this sort is 

G -»• H x U(l) -»• H x Z N . (1) 

For a semi-simple group G, the first of these phase transitions gives rise to monopoles, 
and at the second phase transition each monopole gets attached to N strings. For N > 3, 
a single infinite network is formed which permeates the entire universe ||. The magnetic 
fluxes of monopoles in the network are channelled into the strings that connect them. But 
the monopoles typically have additional, unconfined magnetic charges. The gauge fields 
associated with these charges may include electromagnetic or gluon fields, but may also 
correspond to broken gauge symmetries and have non-zero masses. 

The cosmological evolution of monopole-string networks has been discussed in Ref. 
where it was argued that the networks evolve in a self-similar manner, with a charac- 
teristic scale growing proportionally to cosmic time t. The typical energy of the monopoles 
also grows like t, and the monopoles become highly relativistic. Radiation of gauge quanta 
by accelerating monopoles is the main energy loss mechanism of the networks. Relativistic 
particles emitted by the monopoles (or the decay products of these particles) contribute 
to the spectrum of high-energy cosmic rays. 

For the symmetry breaking sequence 

G^Hx U{1) -> H, (2) 

which can be regarded as a special case of (|I|) with N = 1, the cosmological evolution is 
quite different. In this case, each monopole is attached to a single string, and the second 
phase transition results in monopole-antimonopole (MM) pairs connected by strings. If 
both of the phase transitions in @ occur during the radiation era, then the average 
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monopole separation is always smaller than the Hubble radius, and when MM pairs get 
connected by strings and begin oscillating, they typically dissipate the bulk of their energy 
to friction in less than a Hubble time 

A more interesting possibility arises in the context of inflationary scenario, when 
monopoles are formed during inflation but are not completely inflated away. Strings can 
either be formed later during inflation, or in the post-inflationary epoch. In this case, 
the strings connecting MM pairs can be very long. The correlation length of strings, £, 
can initially be much smaller than the average monopole separation, d; then the strings 
connecting monopoles have Brownian shapes. But in the course of the evolution, £ grows 
faster than d, due to small loop production, and to the damping force acting on the strings. 
Eventually, £ becomes comparable to the monopole separation, and we are left with MM 
pairs connected by more or less straight strings. At later times, the pairs oscillate and 
gradually lose their energy by gravitational radiation and by radiation of light gauge 
bosons (if the monopoles have unconfmed magnetic charges). When the energy of a 
string connecting a pair is dissipated, the monopoles and antimonopoles annihilate into 
relativistic particles. 

Although monopole-antimonopole pairs connected by strings do not typically survive 
until the present time, they can produce a characteristic feature in the spectrum of the 
stochastic gravitational wave background 0. The magnitude of this feature depends in 
particular on the intensity of the gauge boson radiation, which determines the lifetime of 
the pairs. 

The production of high energy particles by topological defects is interesting as a 
mechanism of emission of ultrahigh energy particles, in particular photons, nucleons and 
neutrinos. The energies of such particles can be much higher than have been observed 
until now. Particle production by cosmic strings and by annihilating monopole- 
antimonopole pairs || has been previously discussed in the literatute. Here, we shall 
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calculate the radiation power and the spectrum of gauge bosons, both massless and mas- 
sive, emitted by an oscillating monopole-antimonopole pair. If the magnetic charge of the 
monopoles is of electromagnetic nature, then the radiation can be treated in the same way 
as that of an accelerating charge, and our problem reduces to a problem in classical elec- 
trodynamics. In the case of a chromomagnetic charge, the picture is more complicated. 
The monopoles emit heavy virtual gluons which originate a parton cascade. The hadrons 
are produced in the usual way beyond the confinement radius. The decay of pions results 
in high energy neutrino radiation. The heavy gauge bosons, such as W and Z, can also 
be produced for some range of parameters of the MM pair and the string. 

The problem studied here is in fact a general theoretical problem: the radiation of 
heavy particles by an accelerating monopole. Our results can be reformulated in terms of 
the basic parameters of a shock problem, the proper acceleration of the monopole. 

The dynamics of MM pairs connected by strings is rather complicated and has 
not been studied in any detail in the litterature. We shall therefore concentrate on the 
simplest case of an oscillating pair connected by a straight string, for which the equations 
of motion can be solved exactly. In addition, we shall consider the case of a harmonic 
oscillator motion of the monopoles which has some of the features expected in more general 
monopole-string configurations. Taken together, the analyses of these two cases will allow 
us to draw a reasonably complete picture of the radiation spectrum in the general case. 

After reviewing the general formalism for gauge boson radiation in the next Section, 
we calculate the radiation spectrum of massless gauge bosons in Section [Til]. Radiation 
of massive gauge bosons is analysed in Section [IV]. The spectrum of hadrons resulting 
from radiation of gluons and their subsequent fragmentation are discussed in section |V|. 



Finally, Section [VT| contains some concluding remarks. 
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II. GENERAL FORMALISM 



The characteristic monopole radius 5 m and string thickness 5 S are determined pri- 
marily by the corresponding symmetry breaking scales, r\ m and r\ s . Typically, 5 m ~ 7/" 1 
and 5 S ~ r]' 1 . We shall assume that the two symmetery breaking scales are well sepa- 
rated, 7] s <C r] m , and thus the monopole radius is much smaller than the string thickness, 
5 m < S s . 

Assuming also that the string length is much greater than its thickness, we can treat 
monopoles as point particles and strings as infinitely thin lines. The monopole equations 
of motion can then be written as H 

— = —X v , (3) 
ds 2 m 

where /i is the mass per unit length (and tension) of the Goto-Nambu string, m is the 
monopole mass, r is the proper time of the monopole, and X u is a unit vector orthogonal 
to the monopole wordline and oriented into the string worldsheet. It follows immediately 
from Eq. © that 

a = £, (4) 
m 

is the proper acceleration of the monopoles. Eq. can also be rewritten in a 3- 
dimensional form: 

7 3 x(i) = a 7 (s) n. (5) 

where, 7 = (1 — x 2 ) -1 / 2 is the Lorentz factor of the monopole, 7^ is the Lorentz factor 
of the string at the location of the monopole, and n is a unit vector pointing from the 
monopole in the direction of the string. Eqs and © have a simple physical meaning: 
in the instantaneous rest frame of the monopole, the magnitude of its acceleration is given 
by Eq. (|J) and the direction is given by the direction of the string. The motion of the 



string is described by the usual Goto-Nambu equations. A simple solution of Eq. @ 



describing an oscillating MM pair connected by a straight string is |10 



x(t) = ± 



sgn(t) 



7o - V 1 + (7o v - a\t\) 2 



(6a) 



a 



y(t) 



z(t) = 0, 



(6b) 



where Vq and 70 = (1 — v^)^ 1 ^ 2 are respectively the maximal velocity and the maximal 
Lorentz factors of the monopoles, reached at t — 0. The upper and lower signs correspond 
to M and M respectively. 

The solution @ is valid for \t\ < r y vo/a. At t = —■yoVo/a, the monopoles are at 
rest, with the string having its maximum length, L = 2(70 — l)/a. At t — +70^0/0, 
the monopoles come to rest again, with their positions interchanged. Eq. describes 
only half a period, T = A^qVo/cl, the other half period being obtained by exchanging the 
positions of the monopole and antimonopole. A peculiar feature of this solution (^) is 
that the monopole accelerations abruptly change direction when the monopoles meet and 
pass one another. Between these sign changes the monopoles move at constant proper 
acceleration a. The monopole and antimonopole meet at t — and could be expected 
to annihilate. However, this solution is considered as an approximation for an almost 
straight configuration where the monopole and antimonopole would merely come close to 
each other and would not collide. Besides, as we already mentioned, the monopole radii 
are much smaller than the string thickness, thus the monopoles are not likely to collide, 
even for a straight string. 

To find the radiation power of gauge bosons by accelerating monopoles, we first 
consider radiation by ordinary gauge charges. The gauge field A^(x) produced by the 
4-current j^{x) is determined by the equations 



where F va = d v A a — d a A v is the field strength and M is the gauge boson mass. The 
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F v \ v + M 2 A a = j a , 



(7) 



current corresponding to a pair of equal and opposite point charges is 

f(x) = qx u + 5 {3 \x - x+(t)) - qx"_5 {3) (x - x_(t)), (8) 

where x + (i) and x_(£) are the trajectories of the charges and = dx u ±jdt. In the case 
of a periodic motion, the radiation power can be found from the following equations 

P = Y,Pn = Efd^, (9) 

§--^;K,^K.k). do. 

Here, dP n /dQ is the radiation power at frequency u n = 1-nnjT per unit solid angle in the 
direction of k, |k| = (u; 2 — M 2 ) 1 / 2 , T is the period of the oscillation, and 

1 r T 



n,k) = — J dtexp(iu n t) J d 3 xexp(-ik-x)j fJ -(x, t) (11] 



is the Fourier transform of the current. Eq. ( |I0D is derived in the Appendix The 
derivation is similar to that of the gravitational radiation power in Weinberg's book [ |TT|j . 

Because of electric-magnetic duality, the same equations apply to radiation by mag- 
netic monopoles. Hence, the radiation power from an oscillating MM pair can be found 
from Eqs. (@)-(0) with q interpreted as the magnetic charge. 

For a one dimensional motion, the current has only two non-zero components: j° 
and j 1 . It is also conserved, V„j u = 0, which in Fourier space can be written as 

UJnf = k x j\ (12) 

where k x = k 1 . Combining this with Eq. (110), we have 



^ = ^(i-%iK> k )l 2 - (13) 



The motion of the MM pair in Eq. ([]) has the properties 
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x+ (t) = -x-(t) = -x+(-t) (14a) 
x+(t) = -x+(t + T/2). (14b) 

In such a case, it is easily seen that the system radiates only in odd modes, and that for 
n odd, we have 

8q r T / A 

ji(o; n ,k) = — / x(t)dtcos(u n t)cos[k x x + (t)]. (15) 
1 Jo 

Before we proceed with the calculation of the radiation spectrum, we shall briefly 
comment on the limits of applicability of the classical treatment of the radiation process 
that we adopt in this paper. The problem of radiation by a monopole pulled by a con- 
tracting string is equivalent to that of a monopole moving in a longitudinal magnetic field, 
or of an electric charge moving in a longitudinal electric field. Such processes in strong 
fields are characterized [12|, [13|] by the invariant x given by 



X ={qlm A )^{F^p u f = f/m\ (16) 

where q, m and p u are respectively the charge, mass and momentum of the radiating 
particle and / is the force acting on the particle in its instantaneous rest frame. In the 
case of a monopole connected to a string the corresponding parameter is x = At/m 2 = 



a/m. The case x ^ 1 corresponds to the quasiclassical regime [ 12fl . Hence, our classical 
treatment of radiation is justified provided that /i <C m 2 . 

III. RADIATION OF MASSLESS GAUGE BOSONS 

We begin by calculating the radiation of massless gauge bosons by an oscillating 
MM pair described by Eq. (Q). In this case, Eqs. (|i~3|), (|I~5l) for the radiation power can 
be transformed to 

Pn = A 1 ~ U 2 )[S n (u)] 2 du, (17) 



S n (u) = / o - "° cos«) sin M (_ - ej - __ ) de. (1* 



where we have introduced new integration variables u = k x /u n and £ = niiax + {t)/2'j v . 
These integrals cannot be evaluated analytically. However, the high-frequency asymptotic 
behaviour can be obtained as an expansion in powers of 1/n. This can be done simply by 
repeatedly integrating by parts in Eq. (]T8|) |10[ . 



The leading term of the expansion is 



SnM « ^T^SI (19) 
nn% (1 — u 1 VqY 



and the corresponding power spectrum is 



where 



P^A(m\ (20) 



A - — ['(I - x2) [7o 2 + 3* 2 (7o 2 -l)] 2 ix f2n 



For ultrarelativistic monopoles, 7 3> 1, this can be approximated as 

(22) 

The 1/n expansion is valid as long as the first neglected term is small compared to the 



57T 3 V n ) 



leading term (^D- A somewhat lengthy but staightforward calculation gives the condition 
n ^> 7q. At lower frequencies, the n-dependence is rather complicated and analytic 
approximations are difficult to obtain. 

We computed the radiation spectrum numerically by integrating Eqs. Ql7|) , (fl8|) 
for various values of 70. A generic example of this spectrum, obtained for 70 = 10, is 
plotted in Figure p]. It exhibits the expected behaviour (p2] ) at high frequencies and is 
approximately flat at lower frequencies. The transition happens around n ~ 7q which is 
the characteristic frequency ZJ ~ 7 2 cui, at which most of the power is radiated. We note 
that in the rest frame of the monopole it is 

u res t ~ cU/70 ~ 7o^i ~ a. (23) 



The total power radiated by the system can be evaluated directly using the standard 
formulas of electrodynamics, 

P= { ^f du(l - it 2 ) f ^([ 7o V - ^oO 2 + 1 - u 2 }~ 3/2 + [7oV + ^oO 2 + 1 - «T 3/2 ) 2 - 
An Jo Jo 

(24) 

When the monopole and antimonopole reach ultrarelativistic speeds, their radiation be- 
comes focused in their respective forward directions, which are opposite, and thus inter- 
ference between the monopole and antimonopole fields is negligible. This means that the 
total power radiated by the system is close to twice the power radiated by the monopole 
alone. The radiation intensity of a single charge q moving with proper acceleration a is 
TA\ P = q 2 a 2 /Gtt. Thus, for high velocities of the monopole and antimonopole 7 3> 1, 



the total power tends to 

P^i-(ga) 2 . (25) 

The numerically calculated power has been plotted as a function of 70 in Figure || and 
indeed quickly tends to this asymptotic limit. 

The monopole trajectory (|6]) that we studied so far is rather special. As we already 
mentioned, the monopole accelerations change discontinuously as the monopoles pass one 
another. It is this feature that is responsible for the power-law (rather than exponential) 
falloff of P n at large n. To see how the character of the spectrum is changed in a more 
generic case when the acceleration varies smoothly, we considered a simple harmonic- 
oscillator-type motion of the monopoles, 

x + {t) = Xq sin fit = —X-(t). (26) 

The period of this motion is T = 2ir/Q, the maximum velocity reached by the monopoles 
is Vq = f2xo, and the maximum Lorentz factor is 70 = (1 — v 2 ) -1 ^ 2 - The parameters Xq 
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and f2 should of course satisfy v = £lx < 1, but we can shoose v to be arbitrarily close 
to 1. The proper acceleration of the monopoles in Eq. (^) is 

_ tt 2 x sin Ot 

Its maximum value determined by condition a(t m ) = is reached at sinf2t m = 
±1/V2j v , if V > l/y/3. It is given by 

2 =7o^ (28) 



m 3v/3 

For ultrarelativistic monopoles, 7 ^> 1, the acceleration remains of this order of 
magnitude, a ~ 75^, only for a brief interval of time, 

/ \ 1/2 

At- ^ -(loty- 1 . (29) 
\\a m \J 

Since P oc this part of the oscillation period gives the dominant contribution to the 
radiation power in the massless case. 

The power radiated in massless gauge bosons by the harmonic solution (|26|) can be 
computed from the standard formulas ( |13"D and (|15|) . It gives 



P n = -{nnqf f 1 du{— 2 - l)J 2 n {nuv ), (30) 

7T JO U 



where J n is the n-th Bessel function. 

For large frequency modes n>l and a parameter x < 1, the Bessel function can 
be expanded using the formula 

1 e -2n(l-x 2 ) 3 / 2 /3 

J 2 Jnx)~ . (31) 

to get the simplified expression for the power 

P B «np) / d«(^-l) 6 , ■ (32) 
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At high frequencies, the radiation is focused in the monopoles forward direction, which 
corresponds to u ~ 1. This translates in (|32"D to the steep fall off of the exponential. 

For not too high frequency modes 1 C n ^ 7q, the expression of the Fourier 
transform of the current can be somewhat simplified 

8q f T / A 

j (nQ,k) — — dt sin(nfii) sin[nu(l — I/70) sin(f)t)], 
T Jo 

~ — / dvsininv) sin [rut sin i? , (33) 



7T JO 

where u = k x /Q. This corresponds basically to taking v o = 1 in j - Thus, in this frequency 
range, the approximation fl3"2|) takes the form 

P„ „ „ * VT^exp (-|(1 - . (34) 

A simple change of variable £ = 2n(l — w 2 ) 3 / 2 /3 enables us to evaluate the integral: 

(fig) 2 /"W3 ^e"« 



2tt 2 7o [1 - (3£/2n) 2 / 3 ] 3 / 2 

WrV^W (35) 



2tt 2 7o 2tt 2 

So the spectrum is expected to be flat in this frequency range. 

For very high frequency modes n ^> 7$, we can expand the argument of the expo- 
nential in (^) around u — 1, and simply take all the other slower varying terms as being 
constant. This gives 

P n « ^(fig) 2 ne- 2n/3 ^ 3 J + °° vdve~ 2nv > ho 

(36) 

27r 2 t>o w 

The radiation spectrum (^) was computed numerically for various values of 70. A 
generic example of this spectrum obtained for 70 = 25 is given in Figure |3|. As found in 
(|3~5l), the low frequency behaviour is flat, up to frequencies uj ~ a m 7o- The expected high 
frequency behaviour ( p6[) then sets in at higher frequencies. Compared to the spectrum 
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obtained for the exact solution @, this spectrum is very similar: flat at lower frequencies 
then decaying fast. But, for the harmonic solution, the flat part has a smaller value, 
extends farther in frequencies and the fall off at higher frequencies is much faster. 

The total power radiated by the two monopoles is easy to evaluate when they are 
ultrarelativistic 70 ^> 1. In this limit, as explained previously, interference between the 
monopole and antimonopole fields is negligible and the total power is just the sum of the 
power emitted by the monopole and antimonopole separately. The radiation intensity for 



a single charge moving with a proper acceleration (E71) is [14] P = 9q a^/1670. Thus, for 



ultrarelativistic monopoles 70 3> 1, the total power tends to 

P - = kW- (37) 

870 6 

Contrary to the previous case, Eq. (125|) , this has an extra factor 70 in the denominator. 
It comes from the fact that for ultrarelativistic monopoles, the main contribution to the 
total power comes from the interval of time At in (|29|) when their acceleration is near its 
maximum a m so that P ~ q 2 a 2 n At/T ~ q 2 a 2 n / r )Q. This total power is also consistent with 
the picture of a spectrum consisting of a low frequency flat behaviour P n ~ A^q A up to 
n ~ 7g as shown in (^), and then a negligible exponentially decreasing high frequency 
behaviour which can be ignored. This gives a useful approximation for the spectrum of 
massless gauge bosons, 



(Qa) 2 n 2 

P(E) = ^-^ = P for E < - 7o 3 fi = E , (38a) 

2tt z 3 

P(E) = for E > E , (38b) 



which has the same total power as the exact spectrum and will be a good approximation 
for large 7-factors. 

It is instructive to compare the results we obtained in this Section with the case 
of synchrotron radiation by a charge moving at a constant speed in a circular orbit [ [14}| . 
The power spectrum in this case exhibits a peak around u ~ 7 3 ^, where 7 is the Lorentz 
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factor and f2 is the angular frequency of the rotating charge. The proper acceleration of 
the charge is a w j 2 fl, so once again the characteristic frequency is uj ~ 7a. 



IV. RADIATION OF MASSIVE GAUGE BOSONS 

We will now consider the radiation of massive gauge bosons of mass M. For sim- 
plicity, we will consider only the harmonic oscillation solution ( p6|) in this section. 

The power radiated in massive gauge bosons is still given by the general formula 
(O) but now k = u n u n , where we have introduced 



u n =y/l-AP/u^ (39) 

and uj n = nQ. The current and its Fourier transform remain the same as in the massless 
case, and thus we obtain a very similar formula 

P n = r du{- 2 - l)J 2 n (nuv ), (40) 

TT JO U 

where only the upper bound in the integral has changed. The massless gauge boson 
radiation (f4~0"D is obviously recovered in the limit M — > or n — > 00. Using the asymptotic 
expression of the Bessel function ([H]), valid for n 3> 1, we obtain 

P. M f *L_ ( Izjd) exp - AS)^) (41) 

In the case of massive gauge bosons with M 3> a m , where a m is the maximal acceleration, 
both functions in front of the exponent can be taken at value u = u n and after integration 
we obtain 



where 
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The function (fi3h has a minimum 



/-VSg (44, 



which is reached at 



n = n m = v^-^ — • (45) 
The second derivative of /(n) at n = n m is 

2 n 

fnn = T^it^M' (46) 

The spectrum ( f4~2"D has a maximum at n ~ n m . 

At very large n S> (M/a m )n m we obtain the spectrum (37) of a massless caseQ as 
it must be according to (|40|). 

At n rn <C n < (M / a m )n m the energy flux is exponentially suppressed as 

q 2 7o 2 M 2 1,2 



At M/Q Cn< n m the energy flux is also exponentially suppressed as 

Pn **&% exp(-^M^ 3 /^ 3 ). (48) 
In the region of the maximum n m , the spectrum is 



q 2 I / -v M\ (l f( ,2 

12^ 6XP {'^W J ^ \2 fnnin " ^ 



(49) 



Therefore, in the case we are considering M ^> a m , the characteristic energy of emitted 
gauge bosons is 



1 This can not be seen directly from (|42| ) as the dominant term in this regime was already 
neglected. To get this extra term, the factor (1 — u 2 ) in (|4l|) must be kept and integrated with 
the exponential term 
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E 7H n m Q « 7o M, (50) 

while in the massless case it is E w 7oa m . 

The power spectra ([40"D have been computed numerically and are shown in Figure |] 
for various values of M/a m . Note that our analytical formulae given above are valid only 
for M/a m ^> 1. One can see that at large values of n these spectra tend to the massless 
spectrum shown by a solid line. 

The total radiated power, obtained by integrating ( fiEf ) 

/2tt q 2 7^ 2 r—— — / 2Mv \ 

is exponentially suppressed in the case M 3> a m . 

The formulae above are valid for harmonically oscillating monopoles with a fre- 
quency Q. We can reformulate these results for the instantaneous radiation power from 
a monopole moving with a proper acceleration a and Lorentz factor 7. This is possible 
to do because for M 3> a m the radiation is strongly dominated by four very short bursts 
during the period T. The time dependence of the radiated power is determined mainly 
by the factor 

( 2 Mv\ 

p(t)ocexp \-?>W))' 

and the bursts occur at the maxima of the proper acceleration a(t). Expanding a(t) near 
the maximum, a(t) ~ a m + \a m {t — t m ) 2 , and using Eq . ([27D we obtain an estimate for the 
duration of the bursts 

^ / a^^A 1/2 ^ _J_ /a^N i/> 



It is (a m / Mvq) 1 / 2 times shorter than the corresponding duration in the massless case (|29|) . 

The fact that the dominant part of the radiation comes out in short bursts at a m a m 
suggests the following ansatz for the instantaneous power: 

P(a m , 7 ,£) « k^P har (E,n). (53) 
16 



Here, 7 = y §7m is the Lorentz factor of a harmonically oscillating monopole at the 
moment when it reaches a = a m , k ~ 1 is a numerical coefficient, and Ph ar {E, Q) is given 
by ( |49|) . We expect Eq.([53"D to be valid for E near the maximum of the spectrum at least 
for \E — E m \/E m <C 1. Integrating the instantaneous spectrum ( |53"D over the period T, 
with a(t) and 7(2) for the harmonic oscillator motion, and comparing with the averaged 
oscillator spectrum ( |4*9| ) we determine the normalizing constant, k = 2y / 2/(3y / 7r). The 



total power and the spectrum for a relativistic (7 > 1) monopole moving with proper 
acceleration a and Lorentz-factor 7 are thus given by 

V3 2 w f 2M 



^(«) = ^g'aMexp --— , (54) 



P(a,l,E) = -f-=- (^^exp (-\™) exp (-/^(E - E m f) , (55) 

where 

£ m = ^7^. (56) 
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/£S = (57) 



Eqs.(ES) and (|57D indicate that the gauge quanta are emitted with a characteristic energy 



E ~ 7M in a narrow range 

AE/E ~ (a/M) 1/2 . (58) 

In the case of a vector field ^(x) considered above, there is a compensation in 
the energy-momentum tensor between the two terms corresponding to two nonvanishing 
components of the field (e.g A Q and A 3 ). Such a compensation is absent for a scalar field 
cf>(x) with an interaction L(x) = q(f>(x)j(x), where 

j{x) = J cLt5 w (x-x(t)) (59) 
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and t is the proper time along the trajectory of the source. In this similar 
analysis, for M 3> a, gives the following expression for the instantaneous radiation power 
of a monopole moving with acceleration a and Lorentz-factor 7> 1: 

* = ^v-M-if) ■ < eo > 



V. PRODUCTION OF HIGH ENERGY HADRONS. 

We shall consider here the production of high energy hadrons by an accelerating 
monopole through the most relevant process for this purpose: emission of "heavy" off- 
mass-shall gluons. 

The GUT monopole has |15| a chromomagnetic charge, which is screened beyond 
the confinement radius. The monopole-gluon coupling constant is 1/o.qcd- Production 
of off-mass-shall gluons occurs in the same way as the radiation of massive gauge bosons 
(Section IV). The virtualities of gluons Q 2 higher than a 2 are exponentially suppressed 
and thus the typical vituality and energy of the emitted gluon is a 2 and 7a, respectively. A 
virtual gluon decays into two partons, thus starting a parton cascade, which is described 
by the Gribov-Lipatov-Alterelli-Parisi equation [16] . Beyond the confinement radius, this 



cascade is converted into a hadronic jet. This picture is quite similar to hadron production 
in e + e~ annihilation or in deep-inealastic scattering. We expect it to apply as long as 
a > A, where A = 0.234 GeV is the QCD scale. 

The spectrum of hadrons h radiated by a monopole moving with a proper accelera- 
tion a can be expressed with the help of a fragmentation function, Wh{a, E g , x) as 



fimax 

N h (E h )= 9 dE g N g (E g )W h (a,E g ,x)/E 9 , (61) 

J Eh 



where x = Eh/E g and N g (E g ) is the spectrum of heavy gauge bosons calculated in Section 
V. Actually one can use the massless spectrum when considering virtualities less than a 2 . 
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In the rest frame of the monopole, the typical virtuality and energy of the gluon are, 
respectively, a 2 and a, similar to the virtual photon in e + e~-annihilation. We shall use 
the fragmentation function based on the calculations of ref. ]17| (see also more general 



expression in ref. [18|) which depends only on the parameter x: 



W h (x) = — exp , (62) 



x \ a 

where is the normalization constant, taken from the condition 



o 



xW h (x)dx = f h , (63) 



where fh is the fraction of energy transferred to hadrons of type h, xq = J A/ a, and 



1 2tt 

ff ° = W^' (64) 

The value b = (l/127r)(33 — 2f) determines the QCD coupling constant in the one-loop 
approximation as a" 1 = Mn(a 2 /A 2 ), where / is the number of guark flavors unfrozen at 
a given value of a. Numerically, o"o = 0.121v / ^ln 3 ^ 2 (a 2 /A 2 ). Note that the fragmenta- 
tion function (B^) differs significantly from the function introduced by Hill [H] which was 
subsequently used in much of the literature on cosmic rays from topological defects. 
The fragmentation function Wh{x) has a maximum at 

x = x m = x exp(-o- /2). (65) 

The origin of this maximum is related to the effect of coherent radiation of soft gluons by 
a jet of partons JI7||, ||18|| . At small x the size of emitted gluons is larger than the width 



of the parton beam. A jet therefore radiates soft gluons as a single source with a color 
charge equal to the algebraic sum of the color charges of partons in the beam. 



For ultrarelativistic monopoles, the approximation fl38|) for the spectrum of massless 
gauge bosons can be used instead of the exact spectrum. Then, the integration in Eq. 
(|6l|) for the hadron spectrum can be performed and we obtain, for Eh < Eq, 
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N h {E h ) = — / fc - * , 66 

where we have introduced the auxiliary function 

/•+00 „ 

<j>{x) = / e~* dt. (67) 

We obviously have Nh(Eo) = 0, while the low energy asymptotic behavior of the spectrum 
is: 

P 0f 0(l n(xo )/^) Eo 

N h {E) « / — . 68 

a^o 0(V^o/ 2 ) ^ 

A numerical example for the hadron (nucleon) spectrum, which corresponds to a = 
100 GeV and a = 10 TeV^ is shown in figure |5|. The value of fh can be taken as 0.1, 
which is roughly valid for nucleons. The shape of the spectrum depends only on the 
parameter a. 



VI. CONCLUDING REMARKS 

We have analyzed the radiation of massless and massive gauge bosons by accelerating 
monopoles. In the massless case, we calculated the radiation spectrum for an oscillating 
MM pair connected by a straight string and for a monopole undergoing a harmonic 
oscillator motion. Combining this with the standard analysis of the synchrotron radiation, 
we arrive at the following qualitative picture. 

For a monopole of magnetic charge q moving with a proper acceleration a, the total 
radiated power is P = q 2 a 2 /Qn. In our case, a = fi/m, where fi is the string tension and m 
is the monopole mass. In the instantaneous rest frame of the monopole, the characteristic 
energy of the emitted quanta is E rest ~ a, and thus in the observer's frame it is 

E ~ 7a, (69) 
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where 7 is the Lorentz factor of the monopole. For E ^> E the spectral power is strongly 
suppressed, while for E < E the spectrum is somewhat dependent on the type of motion 
(for example, the spectrum is flat for a harmonic oscillator and has a maximum at E ~ E 
for a circular motion). 

The same picture applies in the case of massive gauge bosons, as long as their mass 
is sufficiently small, M < o. In the opposite limit of large mass, M ^> a, the radiation 
power is exponentially suppressed [see Eq. ([)]])], and gauge quanta are emitted with a 
characteristic energy E ~ 7M in a narrow range AE/E ~ (a/M) 1/2 . The interesting 
case for practical applications is the production of hadrons by an accelerating monopole. 
It effectively occurs through radiation of gluons due to the chromodynamic charge of the 
monopole. As in the general case of heavy quanta, the radiation of gluons with virtualities 
\Q\ 2 > a is exponentially suppressed and therefore most gluons are produced with \Q\ 2 < 
a, where the massless limit is approximately valid. Beyond the confinement radius, high- 
energy gluons fragment into hadrons; the corresponding fragmentation function is given 
by Eq.(H). 

The obtained formulae for radiation of accelerating monopole can be straightfor- 
wadly generalized to any other accelerating particle. 

In a monopole-string network, with N strings attached to each monopole, the proper 
acceleration of a monopole is determined by the vector sum of the tension forces exerted 
by the strings. By order of magnitude it is still given by a ~ fi/m. The cosmo logical 
evolution of monopole-string networks is expected to be characterized by a single length 
scale, 

at) ~ Kt, (70) 

where k = const and t is the cosmic time. This scale gives the average distance between 
monopoles and the average length of string segments. The typical energy of a monopole 
is 



E m ~ fj£ ~ fWCt. ( 71 ) 

The corresponding Lorentz factor is 7 ~ (/i/m)/tt, and from Eq.(^) the typical energy of 
quanta emitted at time t is 

E ~ (fi 2 /m 2 )Kt. (72) 

Assuming that radiation of gauge quanta is the dominant energy loss mechanism of 
the networks, one finds that the parameter k is given by 

k ~ yu/am 2 , (73) 

where a = e 2 /47r and e is the gauge coupling. Networks can also lose energy by production 
of closed loops of string and of small nets. The effect of these mechanisms is hard to 
estimate without numerical simulations of the network dynamics. For the time being, 
Eq. (|73f ) should be regarded as a lower bound on k, the upper bound being due to causality, 
£ ^ t. Hence, 

jj/am 2 < k < 1. (74) 

It is clear from Eqs.(|7l|),((|72"D that the energies of monopoles themselves and of the 
quanta they emit can get arbitrarily large at sufficiently late times. In particular, the 
emitted particles can have super-Planckian energies, E ^> m p , where m v is the Planck 
mass. To give a quantitative estimate, we note that the string tension /i and the monopole 
mass m can be expressed in terms of the corresponding symmetry breaking scales as 
fi ~ T] 2 , m ~ 47ir] m /e. These scales are expected to satisfy 

Vew <Vs<Vm< m pi ( 75 ) 

where t]ew ~ 10 2 GeV is the electroweak symmetry breaking scale. It follows from 
Eqs.([F2|),([F4]) that at the present time 
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E/m p > 10 56 (^ 6 /^mJ), (76) 

and thus the emitted particles are super-Planckian provided that 

4N 2 m m p » 1(T 28 . (77) 

This covers a wide range of parameter values. For example, one could have GUT-scale 
monopoles, r\ m ~ 10 16 GeV, and intermediate-scale strings, r/ s > 10 9 GeV, or electroweak- 
scale strings and light monopoles with r\ m < 10 7 GeV . 

The interaction of super-Planckian particles with the microwave background and 
with cosmic magnetic fields, the resulting cosmic ray fluxes, and some other defect mod- 



els that can give rise to such particles will be discussed elsewhere [19]. Conventional 
astrophysical acceleration mechanisms are not capable to account for particles of such 
tremendous energy, and if they are ever observed, it appears that topological defects 
would be the only possible explanation. 



ACKNOWLEDGEMENTS 

One of the authors (V.B.) is grateful to B.R.Webber for providing the manuscript 
of the book [15] before the publication and to B.R.Webber and V.A.Khose for useful 
discussions. The two other authors (A.V. and X.M.) were supported by the National 
Science Foundation. 



APPENDIX A: GAUGE BOSON RADIATION SPECTRUM 

The electromagnetic radiation of a periodic electromagnetic source has been exten- 
sively studied and is well understood ||14j| . However, for the purpose of this work we will 
need a generalization of the standard formulae to the case of radiation of massive gauge 
bosons. 

23 



We consider a gauge boson field of mass M coupled to a source term j^, 

C = - l -F^F^ v + M 2 A* x A[ l + A^f. (Al) 
To this Lagrangian correspond the equations of motion 

(□ + M 2 )A^ = f, (A2) 
and the energy momentum tensor for the bosons 

= -2(d,A x )(d u A x ) - CgT- (A3) 
The power radiated in bosons can be found from 

P{t) = r 2 J dttS r {t), (A4) 

where 

S r = -2(d A x )(d r A x ) (A5) 

is the radial outgoing energy flux and d r is a derivative with respect to the radial distance 
r = |x|. For a periodic motion, this power can be averaged over a period to get 

P = r 2 \% [ cinS r (t), (A6) 



T 

In the case of a massless gauge bosons M = 0, the equation flA^) can be solved 
exactly using the Lienhart-Wiechert potentials, which greatly simplifies the problem. For 
a massive gauge boson there is no such solution. However, Eq. (A"2) is linear and can be 



solved in Fourier space. We therefore introduce the Fourier transforms 



A"(uj n ,k) = ^JdtJ rf 3 xe i( ^*- k ' x) ^(t,x), (A7) 
"K,k) = ^JdtJ rf 3 xe l ^*- k - x )j^(t,x), (A8) 



J 

where T is the period of the source and uj n = 2nn/T. Eq. ( |A2|) becomes 
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(k 2 -u 2 n + M 2 )#Kk) = f(u n , k). (A9) 
This equation is trivially solved to get 

A"(x) = J2 e ~ iU3ntA n(*) ( A1 °) 
n 

An{ * } ~ J (2vr)3 e jfe2-Jfeg±fe' (AU) 



where we have introduced 



k 2 = u£ - M . (A12) 



The expression ([All ) for A^(x) can be greatly simplified in the particular case we are 



interested in, where it is evaluated at a point x far from the source so that at any point 
of the source x', |x| 3> |x'|: 

A,K,x) = J j J d 3 x>r(t,x>) J k2 _ P fc f ± i£ J sinMWe ifc|x - x ' |cose (A13) 

-H irfM « —3^ n ,k n), (A14) 

47r|x — x'| 47rr 

where we have introduced r = |x| and n = x/r. The radial energy flux ( |A5|) can then be 
expressed to first order in 1/r as 

S r = £M*)e <(m ^'j>«, k n)j;(u m , k n). (A15) 

Then the total power ( p\"6| ) becomes 

P = E P - (A16) 



p « = ty/Ki k n)j*(uj n , k n) (A17) 



where P n is obviously the power radiated by the mode n. 
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FIG. 1. A log-log plot of the gauge quanta radiation spectrum P n /(qa) 2 in the case 70 = 10 
(solid line) with its high frequency approximation ( |22| ) (dashed line). The corresponding spec- 
trum for the harmonic solution has been added (dotted line) for comparison. 

FIG. 2. Total power of massless gauge quanta emission Pj (qa) 2 as a function of 70. It quickly 
tends to a limit of 1/37T. 

FIG. 3. A log-log plot of the gauge quanta radiation spectrum P n /(qa) 2 from the harmonic 
solution (^) in the case 70 = 25 (solid line) with its high frequency ( |36[ ) (long dashed line) and 
low frequency (|3^) (dashed line) approximations. 

FIG. 4. A log-log plot of the massive gauge quanta radiation spectrum P n /(qa) 2 for 70 = 25 
and M/a m taking the values 0.1 (long dash), 0.5 (dash), 1 (dash-dot) and 5 (dot). The massless 
case from figure @ has also been added (solid line). 

FIG. 5. A log-log plot of the hadron spectrum (|6^) for a = 100 GeV (solid line) and 
a = 10 TeV (dashed line) with their low energy approximations ( |68[ ) (respectively dotted line, 
dash-dotted line). 
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